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Abstract 

We propose quantum dense coding protocol for optical images. This protocol extends the 
earlier proposed dense coding scheme for continuous variables [S. L. Braunstein and H. J. Kimble, 
Phys. Rev. A 61, 042302 (2000)] to an essentially multimode in space and time optical quantum 
communication channel. This new scheme allows, in particular, for parallel dense coding of non- 
stationary optical images. Similar to some other quantum dense coding protocols, our scheme 
exploits the possibility of sending a classical message through only one of the two entangled 
spatially-multimode beams, using the other one as a reference system. We evaluate the Shannon 
mutual information for our protocol and find that it is superior to the standard quantum limit. 
Finally, we show how to optimize the performance of our scheme as a function of the spatio- 
temporal parameters of the multimode entangled light and of the input images. 

PACS numbers: 03.67.-a, 03.65.Bz, 42.50.Dv 
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I. INTRODUCTION 



The fundamental properties and possible applications of non-classical light have been 
extensively investigated starting from the mid- 70s. In the past decade such novel fields 
of application of non-classical light have emerged, as quantum information and 
quantum imaging . Quantum imaging uses spatially multimode non-classical states 
of light with quantum fluctuations suppressed not only in time, but also in space. The 
promising idea is to introduce optical parallelism inherent in quantum imaging into var- 
ious protocols of quantum information, such as quantum teleportation, quantum dense 
coding, quantum cryptography etc., thus increasing their information possibilities. The 
continuous variables quantum teleportation protocol proposed in 0,0] and experimentally 
realized in [3, 0] has been recently extended for teleportation of optical images in 0, lloT |. 

Quantum dense coding has been firstly proposed and experimentally realized for dis- 
crete variables, qubits, U,Q| and later generalized for continuous variables in Q]. In 

this paper we propose the continuous variables quantum dense coding protocol for optical 

I I 

images. Our scheme extends the protocol [13] to the essentially multimode in space and 
time optical communication channel. This generalization exploits the inherent parallelism 
of optical communication and allows for simultaneous parallel dense coding of an input 
image with N elements. In the case of a single spatial mode considered in 13} one has 
N =1. 

We calculate the Shannon mutual information for a stream of classical input images 
in coherent state. In this paper we assume arbitrarily large transverse dimensions of 
propagating light beams and the unlimited spatial resolution of photodetection scheme. 
That is, we actually find an upper bound on the spatio-temporal density of the infor- 
mation stream in bits per cm 2 ■ s. This density depends on the degree of squeezing and 
entanglement in non-classical illuminating light. Two sets of spatio-temporal parameters 
play an important role in our protocol: i) the coherence length and the coherence time of 
spatially-multimode squeezing and entanglement, and ii) the spatio-temporal parameters 
of the stream of input images. In our analysis we assume that the sender (Alice) produces 
a uniform ensemble of images with Gaussian statistics, characterized by certain resolution 
in space and time (the Alice's grain). 
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We demonstrate that the essentially multimode quantum communication channel pro- 
vides much higher channel capacity than a single-mode quantum channel due to its in- 
trinsic parallel nature. The density of the information stream is in particular limited by 
diffraction. We find that the role of diffraction can be partially compensated compen- 
sation by lenses properly inserted in the scheme. An important difference between the 
classical communication channel (i. e. with vacuum fluctuations at the input of the scheme 
instead of multimode entangled light) and its quantum counterpart is that in quantum 
case there exists an optimum spatial density of the signal image elements, which should 
be matched with the spatial frequency band of entanglement. 

In Sec. II we describe in detail the scheme of the dense coding protocol for optical 
images. The channel capacity of our communication scheme is evaluated in Sec. III. We 
make our conclusions in Sec. IV. The relevant to our analysis properties of spatially- 
multimode squeezing are given in Appendix. 



II. SPATIALLY-MULTIMODE QUANTUM DENSE CODING CHANNEL 



The optical scheme implementing the protocol is shown in Fig. 1. Compared to the 
generic continuous variables dense coding scheme 3] , here the light fields are assumed to 
be spatially-multimode. At the input, the spatially-multimode squeezed light beams with 
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FIG. 1: Optical scheme for spatially-multimode dense coding. 

the slow field amplitudes Si(p,t) and ^(p, t) in the Heisenberg representation, are mixed 
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at the symmetrical beamsplitter BS±. For properly chosen orientation of the squeezing 
ellipses of the input fields the scattered fields Ei(p,t) and E 2 (p,t) are in the entangled 
state with correlated field quadrature components, as illustrated in Fig. 1. 

The classical signal image field A(p, t) is created by Alice in the first beam by means, 
e. g., of the controlled (with given resolution in space-time) mixing device Mod with 
almost perfect transmission for the non-classical field Ei(p,t). The receiver (Bob) detects 
the entangled state of two beams by means of optical mixing on the symmetrical output 
beamsplitter BS 2 and the homodyne detection of quadrature components of the output 
fields Bi(p,t) and B 2 (p, t). This allows for measurement of both quadrature components 
of the image field with effective quantum noise reduction. 

One can give a more straightforward explanation of the sub-shot-noise detection of 
the signal in the scheme shown in Fig. 1. For the symmetrical scattering matrix of the 
beamsplitters 



{Rnm} 



1/1 ^ 



(2.1) 



/ 
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and equal optical paths of two beams, the effective Mach-Zehnder interferometer directs 
the input squeezed field Si(p, t) onto the detector D ± , and similar for S 2 (p, t), thus allowing 
for sub-shot-noise detection of the squeezed quadrature components in both beams. 
The fields at the inputs of the homodyne detectors D x and D 2 are 

B n (p,t) = S n (p,t) + ^A(p,t), (2.2) 

where n — 1,2. In the paraxial approximation, the slow amplitude of light field B n (p,t) 
is related to the creation and annihilation operators b^q, Q) and b n (q,Q) for the plane 
waves with the transverse component of the wave vector q and frequency Q by 

B n (p, t) = -±= £ b n (q, n)e^). (2.3) 
v L 2 T f>n 

In the case of large quantization volume with the transverse and longitudinal dimensions 
L and cT, the summation is performed over the following values of q and Q: q = (q x , q y ), 
q x = -^ n x, Q y = ~J^ n y an d ^ = Tj7 n with n x , n y and n taking the values 0, ±1, ±2, . . .. 
The free-field commutation relations are given by 

B n (p,t),Bl(p',t')} = 5 n , n ,8(p-p')5(t-t'), 
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b n (q,n),bl,(q',n') = S n>n r 5 Q)f y. 



(2.4) 



The value of the irradiance (in photons per cm 2 • s) is equal to Bl(p,t)B n (p,t), and the 
number of photons in the field mode (q, ft), localized in the quantization volume L 2 cT, is 
b^q, Q)b n (q, ft). The observed photocurrent densities 



h(p,t) = BolB^t) + B\(p,t)], 
I 2 (p,t) = B -[B 2 (p,t)-Bl(p,t)}, 



(2.5) 



have the following Fourier amplitudes 
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Bn- 



b 2 (q,n)-bi(-q,-n) 



(2.6) 



where Bq (taken as real) and iBq are the local oscillator amplitudes used in the homodyne 
detection (see the discussion in Subsection IIII B|) . Here and in what follows we denote the 
Fourier amplitudes of the fields and the photocurrent densities by the lower-case symbols. 

The squeezing transformation performed by the optical parametric amplifiers (OPAs), 
illuminating the inputs of the scheme, can be written |3[ as follows: 



s n (q, ft) = U n (q, ft)c n (<f, ft) + V n (q, ft)4(-g, -ft) 



(2.7) 



where the coefficients U n (q,Q) and V n (q, ft) depend on the pump-field amplitudes of the 
OPAs, their nonlinear susceptibilities and the phase-matching conditions (see Appendix 
for definitions of the squeezing parameters). The input fields c n (q, ft) of the OPAs are 
assumed to be in vacuum state. 

After some calculation we obtain for the Fourier amplitudes of the photocurrent den- 
sities: 

i n (q, ft) = B {f n (q, ft) + a n (q, ft)} , (2.8) 

where 



fr(q, ft) = \ e riim cos^i(g, ft) + ie~ ri ^ n) sin ^i(<f, ft) 



-i4>l(q,Q) 



Ci((f,ft) + 



[h.c, (q, ft) -(-& -ft)], 



(2.9) 
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and 



cos^ 2 (g, O) + ie r2im sin fo(q, tt) 



-ic/>2(<f,n) 



[h.c.,(q,n)^(-q,-Q)}, 
represent the quantum fluctuations of the fields at both photodetectors, and 



ai(£fi) = ■j=[a(q,Q;)+a*(-$,-Sl) 



(2.10) 



:2.id 



a 2 (q,ti) 



iV2 



a(q,Q) - a*(—q, — fl) 



are the detected by Bob components of the Alice's signal image. Here a(q, Q) is the 
Fourier transform of classical field A(p, t), defined in analogy to 1)2. 



III. CHANNEL CAPACITY 



A. Degrees of freedom of the noise and the image field 

In order to estimate the channel capacity one has to define the degrees of freedom of 
the noise and the signal in our spatially-multimode scheme. 

We shall assume that all elements of the scheme: OPAs non-linear crystals, beamsplit- 
ters, modulator and CCD matrices of detectors, have large transverse dimensions. The 
squeezed light fields are the stationary in time and uniform in the cross-section of the 
beams random variables. That is, all correlation functions of these fields are translation- 
ally invariant in the p, t space. For the observed photocurrent densities this implies that 
any pair of the Fourier amplitudes (j2.°-j) and (J2.10)) for given (q, ft) and (—q, —ft) result 
from squeezing of the input fields c(q, ft) and c(—q, — fl) and therefore is independent of 
any other pair. 

On the other hand, since the observed photocurrent densities are real the Fourier 
amplitudes i n (q,Q) and v n (—q y — fi) are not independent, while 

i n (q,n)=il(-q,-{l). (3.1) 

For this reason we consider as independent random variables only the noise terms in 
Fourier amplitudes i n (q, O) for Q > 0. The real and imaginary parts of the complex 

6 



amplitudes i n (q, fi) for Q > are related to the amplitudes of the real photocurrent noise 
harmonics ~ cos(<f • p — Qt) and ~ sin(g • p — Qt), directly recovered by Bob from his 
measurements. 

The Fourier amplitudes of the photocurrent densities ()2.8|) satisfy the relation (|3.1|) and 
therefore it is sufficient to take into account only Q > 0. The random signal sent by Alice 
is assumed to be stationary and uniform in the cross-section of the beams. The amplitudes 
a n (q,Q) for Q > 0, n = 1,2, are taken as independent complex Gaussian variables with 
variance a A (q,Q) depending on (q,fl). Since the transformation (j2.11|) is unitary, the 
Fourier classical amplitudes a(q, Q) for any (q, Q) are also statistically independent, and 
the quantity 

a A (q,Q) = (\a(q,n)\ 2 ), (3.2) 

is the mean energy of Alice's signal wave (q, Q) in the quantization volume, where 
a A (q,fl) = <r A (—q,—Q). Here the statistical averaging is performed with the Gaussian 
complex weight function 

^ s »^^r p {-^m}- (3 - 3) 

In what follows we assume Gaussian spectral profile of width q^ for the ensemble of 
input images in spatial frequency domain, 

, P ( ql + ql\ [ i/fix M < fix/2, 

7r(g A /2) 2 V (W2) 2 / [ \Q\ > Q A /2, 

and, for the sake of simplicity, the narrow rectangular spectral profile U(Q) of width Q A 
and height I/Qa in the temporal frequency domain. Since 

E^(?,n) =L 2 TP, (3.5) 

q,0, 

the total average density of photon flux in the image field per cm 2 ■ s is P. The variances 
of the observables i n (q, fi) are finally found in the form 

(i {i n (q,n),ii(q,Q)} + ) = B 2 [a» A (q,n)+a A (q,n)} , (3.6) 

where { , } + denotes the anticommutator. The quantum noise variances in both detection 
channels are given by 

*n A (m = (Ufn(mji(m} + ), (3.7) 



af A (q,fl) =e 2ri ^' fi )cos 2 Vi(g,fi) + e~ 2ri ^' fi )sin 2 Vi(g,fi), (3.8) 

a* A (q, n) = e~ 2r ^ ^ cos 2 4> 2 (q, Q) + e 2r ^ ^ sin 2 ^ 2 (g, fi). (3.9) 

Using these results we can evaluate the Shannon mutual information for our dense coding 
scheme. 



B. Shannon mutual information for the spatially- mult imode dense coding chan- 
nel 

It is well known that in the case of single-mode squeezed light field the statistics of its 
quadrature amplitudes are Gaussian and can be characterized, e. g., by a Gaussian weight 
function in the Wigner representation. In the homodyne detection of squeezed light, the 
statistics of the photocounts are also Gaussian due to the linear relation between the field 
amplitude and the photocurrent density. The discussion of the homodyne detection in 
terms of the characteristic function can be found in Q|. Some considerations for the 
homodyne detection of spatially multimode fields are presented in [lol ]. 

In our quantum dense coding scheme the statistically independent degrees of freedom of 
the noise and the signal are labeled by the frequencies (q, Q) for Q > 0. One can consider 
our quantum channel as a collection of the statistically independent parallel Gaussian 
communication channels in the Fourier domain. The mutual information between Alice 
and Bob for given detector and frequencies (q, fl) is defined as 



if (q,Q) = H*(q,n) - H^ A \q,Q) . (3.10) 
Here H B (q, Q) is the entropy of Bob's observable, and 



H^ A \q,Q) 



is the averaged over the ensemble of Alice's signals entropy of noise, introduced by the 
channel Q]. For a Gaussian channel the mutual information is given by 



a 



A i 



I*(q,n) =ln(l+ _ B ^^ ). (3.11) 



(7 



The quantum noise suppression within the frequency range of effective squeezing and 
entanglement increases the signal-to- noise ratio at the right side of ()3.11j) . The total 
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mutual information I s , associated with the large area L 2 and the large observation time 
T, is defined as a sum over all degrees of freedom and is related to the density of the 
information stream J in bits per cm 2 ■ s: 

I S = E I*(m = L 2 TJ, (3.12) 

n,q,n>0 

where 




(3.13) 



For qualitative and numerical analysis it is natural to associate such quantities as the 
density of the information stream and of the photon flux with the physical parameters 
present in our quantum dense coding scheme. Squeezing and entanglement produced by 
the type-I optical parametric amplifiers (OPA's), are characterized by the effective spectral 
widths q c and Q c in the spatial and temporal frequency domain. The coherence area in 
the cross-section of the beams and the coherence time are introduced as S c = (27i/q c ) 2 and 
T c = 27r/Q c . For simplicity we assume that both OPA's have the same coherence area and 
coherence time. The correlation area Sa and the correlation time Ta of non-stationary 
images sent by Alice, are related to the spectral widths of the signal qA and Qa by 
Sa = (2ir/qA) 2 and Ta = 2tx/VLa- We consider the broadband degenerate collinear phase 
matching in the traveling-wave type-I OPA's. The coherence time T c of the spontaneous 
downconversion will be typically short compared to the time duration Ta of the Alice's 
movie frame. 

The dimensionless information stream J and the dimensionless input photon flux V 
are defined by J = S C T A J, V = S C T A P- That is, we relate both quantities to the time 
duration of the Alice's movie frame and the coherence area of squeezing and entanglement. 

The optimum entanglement conditions in the OPA's are given by 

ri(q,Q) = r 2 (q,£l) =r(q,Q), 

M&ty = ^2(g,fi)±7r/2^(g,fi), (3.14) 
^(0,0) = tt/2. 

We find the dimensionless information stream J in the following form: 

J = jm, {l + V aBA \^ Q) j(d ; /2)2 exp (-g-t|) } , (3.15) 
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where 



a 



RA 



X 0) = e 2r ^' °) cos 2 ^(k, 0) + e~ 2r ^ °) sin 2 ^(k, 0), 



(3.16) 



and dimensionless spatial frequency is defined as k = q/q c . The relative spectral width 
of the Alice's signal d& = qA/q c = {S c / Sa) 1 ^ 2 can be interpreted as the number of image 
elements per coherence length, i. e. the relative linear density of image elements. In what 
follows we assume a simple estimate q c /2 - 



'2k/l, related to the diffraction spread of 
parametric downconversion light inside the OPA crystal, where k is the wavenumber and 
I is the crystal length. 

Quantum noise in the dense coding scheme is effectively reduced for optimum phase 
matching of squeezed beams. As shown in ^, 16, 17 1, an important factor is the spatial- 
frequency dispersion of squeezing, that is, the g-dependence of the phase of the squeezed 
quadrature component. This dependence is due to the diffraction inside the OPA. A 
thin lens properly inserted into the light beam can effectively correct the (^-dependent 
orientation of squeezing ellipses, as illustrated in Fig. 2. 




FIG. 2: Spatial-frequency dependence of squeezing ellipses without (a) and with (b) phase 
correction by means of a thin lens. The degree of squeezing is exp[r(0, 0)] = 3. 



The improvement in the signal-to-noise ratio for different spatial frequencies can be 
characterized by the inverse noise variance a BA {R ) 0) shown in Fig. 3. As seen from this 
figure, the phase correction by means of a lens allows for the low-noise signal transmission 
within the spatial-frequency band of the effective squeezing. 

In our plots for the mutual information density J we keep constant the coherence area 
S c , the degree of squeezing r(0, 0), and the density of signal photons flux. The dependence 
of mutual information density on the relative linear density of the image elements is 

10 







1 K 2 



3 



FIG. 3: Inverse noise variance in dependence of the spatial frequency k for vacuum noise at the 
input (1), and for squeezing with exp[r(0, 0)] = 3 without (2) and with (3) phase correction. 

shown in Fig. 4. For C 1 (large image elements, Sa ^> S c ), the mutual information 
density increases linearly with cIa, since this means improvement of spatial resolution in 
the input signal. Multimode quantum entanglement between two channels of the scheme 
results in much higher channel capacity compared to the classical limit (vacuum noise at 
the input of the scheme). 

On the other hand, for d,A ^> 1 (image elements much smaller than the coherence 
length), the effect of entanglement on the channel capacity is washed out and J goes 
down to the vacuum limit. This is due to the fact that in the limit Sa *C S c almost 
all spatial frequencies of the signal are outside the spatial-frequency band of the effective 
noise suppression, and the channel capacity is finally limited by vacuum noise. 

The phase correction of squeezing and entanglement significantly improves the channel 
capacity, since it brings the spatial frequency band of the effective noise suppression to its 
optimum value. It eliminates the destructive effect of the amplified (stretched) quadrature 
of the noise field at the higher spatial frequencies, as seen e. g. from Fig. 4b. 
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FIG. 4: Mutual information density for the vacuum noise at the input of the scheme (1), and 
for squeezing with exp(r(0, 0)) = 10 without (2) and with (3) phase correction. The density of 
signal photons is V = 1 (Fig.4a), V = 10 (Fig. 4b). 

IV. CONCLUSION 

In this paper we have extended the continuous variables dense coding protocol proposed 
in [3| onto the optical images and calculated the spatio-temporal density of the Shannon 
mutual information. Our multimode quantum communication channel provides much 
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higher channel capacity due to its intrinsic parallel nature. We have considered the 
role of diffraction in our protocol and have found how to optimize its performance by 
means of a lens properly inserted into the scheme. We have shown that, by contrast to 
the classical communication channel, there exists an optimum spatial density of image 
elements, matched to the spatial-frequency band of squeezing and entanglement. 
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APPENDIX A: PROPERTIES OF SPATIALLY-MULTIMODE SQUEEZING 

The main results for spatially-multimode squeezing are summarized in 0. The coeffi- 
cients of the squeezing transformation ()2.7|) satisfy the conditions 

\U n (q,n)\ 2 -\V n (q,n)\ 2 = l, (Al) 
U n (q, tt)V n (-q, -O) = U n (-q, -ti)V n {q, fi), 

which are necessary and sufficient for preservation of the free-field commutation relations 
(|2.4j) . The spatial and temporal parameters of squeezed and entangled light fields es- 
sentially depend on the orientation angle ip n (q,Q) of the major axes of the squeezing 
ellipses, 

MQ, ty = ± arg {U n (q, Q)V n (-q, -Q)} , (A2) 
and on the degree of squeezing r n (q, Q), 

e ±r„(*n) = 1^(^)1 ± \V n (q,n)\. (A3) 

The phase of the amplified quadrature components is given by 

<f> n (q,Q) = -± &Tg {U n {q,n)V:{-q,-ty}. (A4) 
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In analogy to the single-mode EPR beams, the multimode EPR beams are created if 
squeezing in both channels is effective, and the squeezing ellipses are oriented in the 
orthogonal directions. 

For the type-I phase-matched traveling-wave OPAs, the coefficients U (q, fl) and V(q, fl) 
are given by 



U(q, fl) = exp [i [(k z (q, fl) - k)l - 5(q, Q)/2\ } 



coshr(g,fi) + |^sinhr(g,fi) 



V(q,fl) = exp {i[(k z (q,ty - k)l - 5(g,Q)/2]} f ^sinhr(g,Q). (A5) 

Here / is the length of the nonlinear crystal, k z (q, fl) is the longitudinal component of the 
wave vector k(q, fl) for the wave with frequency uj + fl and transverse component q. The 
dimensionless mismatch function S(q, fl) is given by 

6(q, fl) = (k z (q, fl) + k z (-q, -fl) - k p )l w (2k - k p )l + k'^lfl 2 - q 2 l/k, (A6) 

where k p is the wave number of the pump wave, k p — 2k = in the degenerate case. We 
have assumed the paraxial approximation. The parameter T(q, fl) is defined as 



r(q,0) = y/g*-6*(q,n)/4, (A7) 

where g is the dimensionless coupling strength of nonlinear interaction, taken as real for 
simplicity. It is proportional to the nonlinear susceptibility, the length of the crystal, and 
the amplitude of the pump field. 
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